Rejoinder to "Statistical Modeling of Spatial Extremes"

Rejoinder to "Statistical Modeling of Spatial Extremes" by A. C. Davison, S. A. Padoan and M. Ribatet [arXiv:1208.3378].Comment: Published in at http://dx.doi.org/10.1214/12-STS376REJ the Statistical Science (http://www.imstat.org/sts/) by the Institute of Mathematical Statistics (http://www.imstat.org

Statistical Modeling of Spatial Extremes

The areal modeling of the extremes of a natural process such as rainfall or temperature is important in environmental statistics; for example, understanding extreme areal rainfall is crucial in flood protection. This article reviews recent progress in the statistical modeling of spatial extremes, starting with sketches of the necessary elements of extreme value statistics and geostatistics. The main types of statistical models thus far proposed, based on latent variables, on copulas and on spatial max-stable processes, are described and then are compared by application to a data set on rainfall in Switzerland. Whereas latent variable modeling allows a better fit to marginal distributions, it fits the joint distributions of extremes poorly, so appropriately-chosen copula or max-stable models seem essential for successful spatial modeling of extremes.Comment: Published in at http://dx.doi.org/10.1214/11-STS376 the Statistical Science (http://www.imstat.org/sts/) by the Institute of Mathematical Statistics (http://www.imstat.org

Usefulness of the Reversible Jump Markov Chain Monte Carlo Model in Regional Flood Frequency Analysis

Regional flood frequency analysis is a convenient way to reduce estimation uncertainty when few data are available at the gauging site. In this work, a model that allows a non-null probability to a regional fixed shape parameter is presented. This methodology is integrated within a Bayesian framework and uses reversible jump techniques. The performance on stochastic data of this new estimator is compared to two other models: a conventional Bayesian analysis and the index flood approach. Results show that the proposed estimator is absolutely suited to regional estimation when only a few data are available at the target site. Moreover, unlike the index flood estimator, target site index flood error estimation seems to have less impact on Bayesian estimators. Some suggestions about configurations of the pooling groups are also presented to increase the performance of each estimator

Multivariate Copula Analysis Toolbox (MvCAT): Describing Dependence and Underlying Uncertainty Using a Bayesian Framework

We present a newly developed Multivariate Copula Analysis Toolbox (MvCAT) which includes a wide range of copula families with different levels of complexity. MvCAT employs a Bayesian framework with a residual-based Gaussian likelihood function for inferring copula parameters and estimating the underlying uncertainties. The contribution of this paper is threefold: (a) providing a Bayesian framework to approximate the predictive uncertainties of fitted copulas, (b) introducing a hybrid-evolution Markov Chain Monte Carlo (MCMC) approach designed for numerical estimation of the posterior distribution of copula parameters, and (c) enabling the community to explore a wide range of copulas and evaluate them relative to the fitting uncertainties. We show that the commonly used local optimization methods for copula parameter estimation often get trapped in local minima. The proposed method, however, addresses this limitation and improves describing the dependence structure. MvCAT also enables evaluation of uncertainties relative to the length of record, which is fundamental to a wide range of applications such as multivariate frequency analysis

Non-stationarity in peaks-over-threshold river flows:a regional random effects model

Under the influence of local- and large-scale climatological processes, extreme river flow events often show long-term trends, seasonality, inter-year variability and other characteristics of temporal non-stationarity. Properly accounting for this non-stationarity is vital for making accurate predictions of future floods. In this paper, a regional model based on the generalised Pareto distribution is developed for peaks-over-threshold river flow data sets when the event sizes are non-stationary. If observations are non-stationary and covariates are available then extreme value (semi-)parametric regression models may be used. Unfortunately the necessary covariates are rarely observed and, if they are, it is often not clear which process, or combination of processes, to include in the model. Within the statistical literature, latent process (or random effects) models are often used in such scenarios. We develop a regional time-varying random effects model which allows identification of temporal non-stationarity in event sizes by pooling information across all sites in a spatially homogeneous region. The proposed model, which is an instance of a Bayesian hierarchical model, can be used to predict both unconditional extreme events such as the m-year maximum, as well as extreme events that condition on being in a given year. The estimated random effects may also tell us about likely candidates for the climatological processes which cause non-stationarity in the flood process. The model is applied to UK flood data from 817 stations spread across 81 hydrometric regions

On the relationship between total ozone and atmospheric dynamics and chemistry at mid-latitudes – Part 1: Statistical models and spatial fingerprints of atmospheric dynamics and chemistry

We use statistical models for mean and extreme values of total column ozone to analyze "fingerprints" of atmospheric dynamics and chemistry on long-term ozone changes at northern and southern mid-latitudes on grid cell basis. At each grid cell, the r-largest order statistics method is used for the analysis of extreme events in low and high total ozone (termed ELOs and EHOs, respectively), and an autoregressive moving average (ARMA) model is used for the corresponding mean value analysis. In order to describe the dynamical and chemical state of the atmosphere, the statistical models include important atmospheric covariates: the solar cycle, the Quasi-Biennial Oscillation (QBO), ozone depleting substances (ODS) in terms of equivalent effective stratospheric chlorine (EESC), the North Atlantic Oscillation (NAO), the Antarctic Oscillation (AAO), the El Niño/Southern Oscillation (ENSO), and aerosol load after the volcanic eruptions of El Chichón and Mt. Pinatubo. The influence of the individual covariates on mean and extreme levels in total column ozone is derived on a grid cell basis. The results show that "fingerprints", i.e., significant influence, of dynamical and chemical features are captured in both the "bulk" and the tails of the statistical distribution of ozone, respectively described by mean values and EHOs/ELOs. While results for the solar cycle, QBO, and EESC are in good agreement with findings of earlier studies, unprecedented spatial fingerprints are retrieved for the dynamical covariates. Column ozone is enhanced over Labrador/Greenland, the North Atlantic sector and over the Norwegian Sea, but is reduced over Europe, Russia and the Eastern United States during the positive NAO phase, and vice-versa during the negative phase. The NAO's southern counterpart, the AAO, strongly influences column ozone at lower southern mid-latitudes, including the southern parts of South America and the Antarctic Peninsula, and the central southern mid-latitudes. Results for both NAO and AAO confirm the importance of atmospheric dynamics for ozone variability and changes from local/regional to global scales

Extreme events in total ozone over Arosa – Part 1: Application of extreme value theory

In this study ideas from extreme value theory are for the first time applied in the field of stratospheric ozone research, because statistical analysis showed that previously used concepts assuming a Gaussian distribution (e.g. fixed deviations from mean values) of total ozone data do not adequately address the structure of the extremes. We show that statistical extreme value methods are appropriate to identify ozone extremes and to describe the tails of the Arosa (Switzerland) total ozone time series. In order to accommodate the seasonal cycle in total ozone, a daily moving threshold was determined and used, with tools from extreme value theory, to analyse the frequency of days with extreme low (termed ELOs) and high (termed EHOs) total ozone at Arosa. The analysis shows that the Generalized Pareto Distribution (GPD) provides an appropriate model for the frequency distribution of total ozone above or below a mathematically well-defined threshold, thus providing a statistical description of ELOs and EHOs. The results show an increase in ELOs and a decrease in EHOs during the last decades. The fitted model represents the tails of the total ozone data set with high accuracy over the entire range (including absolute monthly minima and maxima), and enables a precise computation of the frequency distribution of ozone mini-holes (using constant thresholds). Analyzing the tails instead of a small fraction of days below constant thresholds provides deeper insight into the time series properties. Fingerprints of dynamical (e.g. ENSO, NAO) and chemical features (e.g. strong polar vortex ozone loss), and major volcanic eruptions, can be identified in the observed frequency of extreme events throughout the time series. Overall the new approach to analysis of extremes provides more information on time series properties and variability than previous approaches that use only monthly averages and/or mini-holes and mini-highs

On the relationship between total ozone and atmospheric dynamics and chemistry at mid-latitudes Part 2

We present the first spatial analysis of "fingerprints" of the El Nino/Southern Oscillation (ENSO) and atmospheric aerosol load after major volcanic eruptions (El Chichon and Mt. Pinatubo) in extreme low and high (termed ELOs and EHOs, respectively) and mean values of total ozone for the northern and southern mid-latitudes (defined as the region between 30 degrees and 60 degrees north and south, respectively). Significant influence on ozone extremes was found for the warm ENSO phase in both hemispheres during spring, especially towards low latitudes, indicating the enhanced ozone transport from the tropics to the extra-tropics. Further, the results confirm findings of recent work on the connection between the ENSO phase and the strength and extent of the southern ozone "collar". For the volcanic eruptions the analysis confirms findings of earlier studies for the northern mid-latitudes and gives new insights for the Southern Hemisphere. The results provide evidence that the negative effect of the eruption of El Chichon might be partly compensated by a strong warm ENSO phase in 1982–1983 at southern mid-latitudes. The strong west-east gradient in the coefficient estimates for the Mt. Pinatubo eruption and the analysis of the relationship between the AAO and ENSO phase, the extent and the position of the southern ozone "collar" and the polar vortex structure provide clear evidence for a dynamical "masking" of the volcanic signal at southern mid-latitudes. The paper also analyses the contribution of atmospheric dynamics and chemistry to long-term total ozone changes. Here, quite heterogeneous results have been found on spatial scales. In general the results show that EESC and the 11-yr solar cycle can be identified as major contributors to long-term ozone changes. However, a strong contribution of dynamical features (El Nino/Southern Oscillation (ENSO), North Atlantic Oscillation (NAO), Antarctic Oscillation (AAO), Quasi-Biennial Oscillation (QBO)) to ozone variability and trends is found at a regional level. For the QBO (at 30 and 50 hPa), strong influence on total ozone variability and trends is found over large parts of the northern and southern mid-latitudes, especially towards equatorial latitudes. Strong influence of ENSO is found over the Northern and Southern Pacific, Central Europe and central southern mid-latitudes. For the NAO, strong influence on column ozone is found over Labrador/Greenland, the Eastern United States, the Euro-Atlantic Sector, and Central Europe. For the NAO's southern counterpart, the AAO, strong influence on ozone variability and long-term changes is found at lower southern mid-latitudes, including the southern parts of South America and the Antarctic Peninsula, and central southern mid-latitudes

Composite likelihood methods for large Bayesian VARs with stochastic volatility

© 2020 John Wiley & Sons, Ltd. Adding multivariate stochastic volatility of a flexible form to large vector autoregressions (VARs) involving over 100 variables has proved challenging owing to computational considerations and overparametrization concerns. The existing literature works with either homoskedastic models or smaller models with restrictive forms for the stochastic volatility. In this paper, we develop composite likelihood methods for large VARs with multivariate stochastic volatility. These involve estimating large numbers of parsimonious models and then taking a weighted average across these models. We discuss various schemes for choosing the weights. In our empirical work involving VARs of up to 196 variables, we show that composite likelihood methods forecast much better than the most popular large VAR approach, which is computationally practical in very high dimensions: the homoskedastic VAR with Minnesota prior. We also compare our methods to various popular approaches that allow for stochastic volatility using medium and small VARs involving up to 20 variables. We find our methods to forecast appreciably better than these as well

A regional Bayesian POT model for flood frequency analysis

Flood frequency analysis is usually based on the fitting of an extreme value distribution to the local streamflow series. However, when the local data series is short, frequency analysis results become unreliable. Regional frequency analysis is a convenient way to reduce the estimation uncertainty. In this work, we propose a regional Bayesian model for short record length sites. This model is less restrictive than the index flood model while preserving the formalism of "homogeneous regions". The performance of the proposed model is assessed on a set of gauging stations in France. The accuracy of quantile estimates as a function of the degree of homogeneity of the pooling group is also analysed. The results indicate that the regional Bayesian model outperforms the index flood model and local estimators. Furthermore, it seems that working with relatively large and homogeneous regions may lead to more accurate results than working with smaller and highly homogeneous regions